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Abstract. The Harper equation describes an electron on a 2D lattice in mag- 
netic field and a particle on a ID lattice in a periodic potential, in general, 
incommensurate with the lattice potential. We find the distribution of the roots 
of Bethe ansatz equations associated with the Harper equation in the limit as 
a = 1/Q — > 0, where a is the commensurability parameter (Q is integer). Using 
the knowledge of this distribution we calculate the higher and lower boundaries 
of the spectrum of the Harper equation for small a. The result is in agreement 
with the semiclassical argument, which can be used for small a. 
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1 Introduction 



The Harper equation (see and also |||7],|8| for recent reviews) 

i> n -i + 2 cos(27ran + 9)ip n + ip n+ i = eijj n , n = . . . , —I, 0, 1, ... , a, 9 E R 

(1) 

in l 2 (Z) describes a quantum particle on the line in two periodic potentials which 
are incommensurate if a is irrational. It is also a model for an electron on a square 
lattice subject to a perpendicular uniform magnetic field (Azbel-Hofstadter). In 
the latter case a is proportional to the value of the magnetic field and 9 is a quasi- 
momentum in the direction along one of the sides of the elementary square (see, 
e.g, jn|). The Harper equation has many deep connections to various fields of 
modern physics, e.g., quasicrystals, quantum Hall effect, and pure mathematics: 
number theory, functional analysis. 

Equation (JJ) is the eigenvalue equation for a particular case of the almost 
Mathieu operator. The spectrum am(a, 9) of the corresponding operator consists 
of Q intervals (bands) if a = P/Q, where P and Q are coprime integers^. The 
spectrum cr^(a) of the Azbel-Hofstadter model is the union over all real 9 of 
0"m(o!,0). At a = P/Q the spectrum a H (P/Q) consists of Q bands. If a is 
irrational, it is known that om(«, 9) is independent of 9, and hence the spectrum 
&h{&) — o"m(«, #)• The plot of cr# (a) as a function of a is called the Hofstadter 
butterfly ||. The case of irrational a is the most interesting. There is a conjecture 
that the measure of the spectrum is zero for all irrational a |9, HID] • This and 



many other questions regarding the spectrum of the almost Mathieu operator at 
irrational a are still unresolved (see |7|, p!2|JT^| for a list of problems). 

Recently a new approach to these problems has been proposed [0,0|. Let 



a = P/Q. First, it can be shown JTJ| that certain Q (one from each band) points 



of the spectrum^ are the solutions e of the following eigenvalue equation: 

i(z~ l + qz)^(qz) - i(z~ l + q~ l z)^> (q^ 1 z) = e*f?(z) (2) 

in the space of polynomials ^(z) = n®=i(. z ~ z k) of degree Q— 1. Here q = e* 7 , 7 = 
ira = nP/Q. This equation is related to the algebra U q {sl-i). One can obtain other 
difference equations similar to (||]) with the same set of eigenvalues. Analogous 
equations are also available for other points of the spectrum ctm(P/Q, 9) fl4] , |T5 



1 That is they do not have a common divisor other than 1. 
2 This is true both for a M (P/Q,6) and for a H (P/Q). 
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In the Appendix we present a simple derivation of (0) which is based on the fact 
that it can be regarded as an equation for the generating function of a finite 
system of orthogonal polynomials related to ([I]). 

Now it is important to note [0 that the following set of relations can be 
obtained from @ by substituting ^(z) = n®=i( z ~ z k) and setting z = z m , 
m — 1, . . . , Q - 1: 

-i + <l gQ Q Y[ qZm ~ Zk , m = l,...,Q-l. (3) 



qz^ + 1 z m - qz k 

Collecting the coefficients at z®~ 1 in (0) gives the expression for the energy: 



-i 



iq Q {q-q- 1 )Y. z m. (4) 



771=1 



Because of the analogy to one-dimensional integrable models, where a similar 
set of equations is obtained, expressions (|3|||) are called Bethe ansatz equations. 
The Bethe ansatz equations for integrable models are solved in the limit of infinite 
system. In our case this limit corresponds to Q — > oo. Thus, we can hope to 
obtain from (|3],|]) information about the spectrum at irrational a. 

The analysis of (|3|,|4]) is only in its beginnings so far. The system ([|,|]) was 
investigated numerically (in the case of e = analytically) in . An analytical 
approach based on the string hypothesis is proposed in |X 



In Section 2 we will obtain a solution (that is the distribution function of the 
roots z^ in the limit as a = 1/Q — > 0, thus explaining some numerical results 
of |1| . To obtain this solution, we will consider besides another set of Bethe 
ansatz equations which follow from (0). 

We find the solution under assumption that for a = 1/Q the roots lie on the 
unit circle, and their limit distribution function p(0) depends piece- wise smoothly 
on the angle 0.0 The function p(0) satisfies an integral equation that follows from 
the Bethe ansatz equations. Besides p(0), another important function, x(0)j 
naturally arises from the Bethe ansatz equations. It differs from unity only when 
p(4>) = const = 1/tt and describes in this case the changing with <ft exponentially 
small in Q correction to the constant distance between neighboring roots. Note 
that a nonconstant p(4>) implies the increase of order 1/Q in this distance with 



3 As usual, Qp{(j))d(j) is the number of roots in the interval 



3 



In Section 3, using the expressions for p(0) and x(0), we obtain the upper 
boundary e max of the spectrum of (1) for small a. The result is e max = 4 — 
2tt/Q + o(1/Q). The value 2ix (—2tc) is the tangent of the lower (upper) boundary 
of the Hofstadter butterfly for small magnetic field. Alternatively, this result 
can be obtained using semiclassical considerations [[17], |18| . Hopefully, the ideas 
introduced in the present work will be helpful in the cases where semiclassical 
techniques cannot be applied. 

2 Distribution of the roots of the Bethe ansatz 
equations for a = 1/Q —> 

In order to obtain the results of this work, we will also need another set of Bethe 
ansatz equations which is derived from (|2|) by first replacing z with qz and then 
substituting z = z m , m = 1, . . . , Q — 1. We get 

Q-i „2 r _ 

II m = l,...,Q-l. (5) 



KQ lz m + 1 2z m) k =l <l Z m ~ Z k 

Henceforth, we will be interested in the simplest case of the Harper equation 
when a = l/Q(P = l, j = ir/Q) and Q is large. Numerical data [|16[] suggests 
that if a = 1/Q all the roots Zk lie on the unit circle, they are simple, and their 
distribution becomes continuous as Q — > oo. Let us therefore assume z\. = e^ k , 
(pk £ R, substitute this expression into and then try to find from these 

equations a piece-wise smooth distribution function p(4>), <f) £ [ — tTj tt] m the limit 
as Q — > oo. After the substitution z k = e l< ^ k , the equations (^|,||) take the form: 

cos(0 m - 7/2) = _ tt 1 sin|Q m - fc + 7) 
cos(0 m + 7/ 2 ) siii|(^ m -0 fc -7)' 

~g = rr 1 Sin |(0 fc - m - 2 7 ) 

2 cos(0 m + 3 7 /2) ^ sin - <f> m - 7) ' U 

We set 4>k+i — 4>k = 8k- According to our assumption 5k is of order 1/Q 
so that p(4>k) = 1/QSk- Smoothness of the distribution function implies that 
(p k — 4>k± P = Tp8k for p = 1, . . . , iV up to at most 0(N 2 /Q 2 ) in the limit as 
N, Q — > 00 in such a way that N/Q — > 0. Let us consider equations (0) in this 
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Figure 1: A few consecutive roots = e lr ^ k of the Bethe ansatz equations near 
the point k = m are shown as they lie with respect to each other for sufficiently 
large Q. The angles S m , 7, A m l are, of course, increased. 

limit. The product is separated into three distinct parts. First, the contribution 
of the two possible points m+n and m+ 2n such that m+n — 4>m — 1 = o(l/Q) 
and </> m +2n — 0m — 27 = o(l/(5), where n is a certain integer which is the same 
in the region of angles 0^ from k = m — N to k = m + N. Define A m j = 
<pm + i - 0m - 7 (See Figure [I]). Then m+2 n - m - 27 = A m+ni „ + A mjT1 . Note 
that if to the main order in 1/Q the distance 5 m 7^ 7, 7/2, 7/3, . . ., then, whatever 
i, both A m>i = 0(1/ Q) and A m+i)i = A mji + o(l/Q). In this case the ratio 
sin |(0 m+2 n - 0m - 27)/ sin \((f) m +n - 4>m - l) is always equal to 2 in the limit 
and, as we shall see, it is not necessary to distinguish it as a separate contribution. 
However, when 5 m = 7,7/2,7/3, ... up to o(l/Q), this ratio is an unknown and 
should be treated separately. Below only the case 8 m = 7 = tt/Q will be relevant 
because in the solution we obtain, 5 m >n/Q for any m . 

The second contribution to the product comes from the rest of the points 
4> m ± p , where p = 1, . . . , N. All the points not included into these two groups give 
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the third contribution. 

Thus, in the limit N,Q — > oo, N/Q — > 0, which we denote just "lim", we can 
write (0) as follows (denoting = lim0 m , e [— 7r,7r]): 

— ^— = lim ABC. (8) 

2cos0 v 1 

Here A stands for the contribution of the points with the indices k = m + n and 
k = m + 2n. Expanding the sinus in l/Q, we can write: 

sin \{n8 m - 27) sin \(<p m +2n - 4>m - 27) 



A 

sin \(4> m +n ~ <Pm - 7) SHI ~(2n5 m - 7) 



(-^ ^m+n,n j ^7 
A m>n / 2n£ m - 7 

The factor S denotes the contribution of the points with the indices k = 
k = m — n, k = m — 2n, and also k = m ± p, p = 1, . . . , N, p ^ n, 2n: 

_ ^ sin \{n5 m + 27) sin \{2n5 m + 27) SL sin \(p5 m - 27) sin \(p5 m + 27) 
sin I (n5 m + 7) sin | {2n5 m + 7) A = A sin | (p5 m - 7) sin | (p5 m + 7) 



l + 2 1 /nd m l-(2 7 /p«J t , 



1 + 7/2n£ m ^ 1 - (7/p^ 



(10) 



Thus 



A m+2n _A 1 - (2 7 /n5 m ) 2 " l-(2 7 /^ 2 



m } 



2 

m I 



(11) 



\ ^ A m+n _! J 1 - ( 7 /2n<5 m ) 2 AA 1 - ( 7 /p$ 

where we used the formula cos(7ra/2) = n&Lo(l — a 2 /{2k + l) 2 ). Finally, C 
in (§) is the contribution of the points with the indices k outside the range 
m — N, . . . ,m + N. We denote the latter restriction by primes at the sums and 
products: 



sin \{(j) m - <p k + 27) TT , 1 + 7 cot l(4> m - (p k ) 

(12) 



AA sin ~(0 m - fc + 7) A fe A 1 + 2cot|(0 m -0& 



exp 



V 2 k 2 
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where the last product was replaced by the exponent of its logarithm. The small- 
parameter expansion of the logarithm is valid here because in this range of k the 
quantity 7cot|(0 m — <pk) is of order 1/N or less. Note that because we take 
the limit N, Q — > oo in such a way that N/Q — > 0, the factor fll2|) includes the 
contribution of all roots of the Bethe ansatz equations except for those lying in 
the infinitesimally small neighborhood of z m = e l< ^ m . Replacing the sum in flT2|) 
by the integral and inserting ([TTD , ( |12|) into (H), we obtain: 



2cos0 \ A m ,„ 

'7r_ r 1 



exp (l^-P- /^ cot ^)P(^)^) • 

A similar analysis of (§) yields the result: 



(13) 



lim 



A m __ n n 



exp (^nV.p. J cot - (0 — cj)p(u;)gL^ . (14) 



We see that the function = lim(A m „/A m _ nin ) is continuous!] Hence, we 
can substitute lim(A m+n n /A m n ) = x(0)> which is expressed from (0), in (PT3|) . 
Thus, 

cosh < — l/.p. / cot -(</> — u)p(u>)du \ cos(7r 2 p(0)). (15) 



4cos0 1 2 7-tt 2 

Substituting here the energy 

e = 2n T e iuj p(tu)dcu, (16) 

we obtain an integral equation for the distribution function p(4>). It is not the 
only possible form of the integral equation for p(0) we can obtain. We can derive 
others by considering other sets of Bethe ansatz equations which follow from (|2[) 
in a similar manner as (]5|). 

Equation fllBT) reflects the well-known symmetry of the roots {-2 m }. They lie 
symmetrically with respect to the real axis. So p(— 0) = p{4>). 

Note that the function p(0) must also satisfy the normalization condition 

p{<f>)d<f>=l. (17) 



We omit the index n of x(</>) because in the solution we will obtain, it is always n = 1. 
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Eliminating the exponents from ([13]) and (|T4| ) we get a simple relation between 
X{4>) and p(4>): 



- e^xW) = 2cos0cos(vr 2 p(0))(l + x(<f>))- (18) 

We can guess the form of the solution p(0) using the following argument. 
Since p(0) is even, it follows from QUI) that x(0) = xC 71 ") — 1- Thus, in some 
neighborhood of the points <f> — 0, tt equation (|18|) simplifies to the form e = 
— 4cos0cos(7r 2 p(0)).0 This can hold only in the range |cos0| > e/4. Outside 
of it, x{4>) can no longer be 1, which means that p(<f>) can be only a piece-wise 
constant function there taking some of the values 1/tt, 2/tt, ... or zero. A precise 
form of p(0) is given by the following 

Statement. p(0) = 4, arccos _ 4 £ COS(t> if 4> E [-if>, if>] U [tt - if), tt] U [-tt, -tt + if>] , 
where ^ = arccos f, < arccos(x) < 7r; e G [0,4]. 
p(0) = 1/tt if <f> G [V, vr/2) U (-tt/2, -V]; 
p(0) = if G [vr/2, tt - if)] U [-tt + up, -tt/2}. (See Figure [|) 

The same but reflected with respect to the imaginary axis distribution corre- 
sponds to the interval e E [—4, 0]. 

Henceforth, we always assume e E [0,4]. 

We verify the statement substituting this expression for p(<f>) in (|15|), fll6|) , 
flTTD and integrating. For flj"6|) we have: 



2tt / e lu} p{u)du = 

J —IT 

2 + / e l "ck> + \ + - e luJ arccos — rfcj 

\ 1/ V J-tt/2) [■/-V' A— J TT — 4cOSCo> , 

A 8 £ J 

4 / cos u arccos duo = 

it Jo 4 cos U) 

2e r u xdx 2ev r°° \fu 2 — v~ 2 du 

4 / arctanx , = = 4 



(19) 



TT JO V^ 2 — X 2 TT Ju- 1 u(l + U 2 ) 

where v = J (4/e) 2 — 1. Henceforth, we take *Jx > if x > 0. 



Comparing this with (|13|) we see that only the points lying in the infinitesimally small 
neighborhood of z m , whose contribution is described by (11), influence p{<f> m ). 



8 




Figure 2: Distribution of the roots of the Bethe ansatz equations in the limit as 
a — 1/Q — > 0. In the upper figure, the regions on the unit circle are shown where 
the roots lie with the density p(4>), where is the angle. In the lower figure, p{<p) 
for some e is plotted as a function of the angle. 
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Figure 3: The contour in the complex plane used to evaluate the integral (p0|). 
The thick lines denote cuts. There are two simple poles at the points /i and — \l. 



Since p(0) is an even function, it is sufficient to verify expression (^) only for 
G [0, 7i"]. There are three distinct cases: G [0,-0], G [i/j,7t/2\, G [7r — i/j,tt]. 
Most of the relevant integrals are easy to evaluate. Perhaps only the calculation 
of the following one is somewhat less straightforward: 



1 r 



ooti( 



uj) + tan -( 



arccos{e/(4cosa;)} , 
V.p. / —- ; duj 

J-if> sm(0 — uj) 

7r 1 + sin 

— m -, 

2 1 — sin d) 



uj, 



i smi 



arccos • 



-dw 



-V.p. 



4 cos a; 
" x arctan x dx 



o x z 



ar 



(20) 



where v = y (4/e) 2 — 1, jti = y (4 cos (ft / e) 2 — 1, cos0 > e/4. The contour used 
for evaluation of the last integral in (|20"D is shown in Figure |3|. 

From (0) we obtain the expression for x{4>) on the set G [V>, 7r/2) U 
(— 7r/2, —-0] where p(0) = 1/tt (note that x(0) = 1 if G [— ijj,ifj] U [7r — -0, 7r] U 
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[-7T, -7T + V]): 



A 1 1 - a/1 - (-cos0) 2 
x(0) = lim-^ = V ; , *e[1>,*/2). (21) 

fl m-l,l 1 + ./l - cos 0)2 



- e 



Thus, in the region £ (i/j,it/2) we have A m -i+fc,i = x(0)( fc / Q ) Q A m _i i i, < 
x(0) < 1, which means that as increases from ip, the distance between the 
roots m approaches the constant value n/Q exponentially fast. We shall use 
the function m the next section to determine the finite Q correction to the 
boundary e = 4 of the spectrum. 

Because of the mentioned symmetry of the Bethe ansatz equations, it is suffi- 
cient to write the expressions for x{4>) only for £ [0, ir]. This symmetry implies 

x(-0) = x(0)- 1 - 

3 The boundaries of the spectrum of the Harper 
equation for small a. 

Consider now the simplest case e = 4 in more detail. In this case p(<fi) = 1/tt if 



£ (— 7r/2,7r/2) and zero otherwise. Equation (21) takes the form 



X(0) = ^-^, 0e[O,vr/2). (22) 
1 + sm <p 

We shall calculate the largest finite Q correction to the energy e — 4. It turns out 
that the shift (from the uniform distribution defined below) of the roots of the 
Bethe ansatz equations described by (j22"D is responsible for this correction. Since 
according to (^3), A m l decreases exponentially as m moves away from the point 
= 0, only an infinitesimally small neighborhood of = gives considerable 
contribution to the shift. Let us assume the following distribution of the roots 
0o = 0, m = mix /Q + YT=o A;,i> m = 1, . . . ,Q/2-l, 0_ m = -0 m (we set Q even 
for definiteness) . The mentioned property of A^j implies that if m — > ^ as 
Q — > oo, then with an exponential accuracy m = mir/Q + (m/\m\)s/Q s , where 
s/Q 5 = J2iZo The quantities s and 5 are to be found from the Bethe ansatz 
equations. Henceforth, we assume m > 0. Further, set 



sin(A m ,i/2) 
sin(A m _ 1>1 /2) 




(23) 
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Retaining now in equations (||) not only the largest in 1/Q quantities (of order 
1) as before, but also the next largest (that is of the order of 1/Q, 1/Q 13 , 1/Q S ), 
we obtain with this accuracy the following expression for the logarithm of (|5|) 

7i Tim irm r m 

- tm —^ xi —) + -- 

2s [ 7r _ 2 ir(m — k) 1 n ^ ' it(m — k) (24) 



Replacing the first sum over k in (|24|) by the integral gives 

7T ° . _ 2 vr(m - jfe) r° . _ 2 - x dx _ 2 

k=-h^) Sm ~ ^ 2 Sm T " l-cos^ + sin^- 

(25) 

As to the second sum over k in fl2"4|) , we need to evaluate it up to 1/Q. Using the 
Poisson summation formula, we have 

7r^/ ^ n(m — /c) 7r QA^ 2m ^ 7r(m — /c) 

^ fc 2Q Q k=-(Q/2-l) 2Q 

/tt i o 7rm | 7r 7T7TT. ,- r , OO / 7 \ 

77 | Q 7Tm | 7T 7T?77- ^ QQ 7T | Q 7TTO j TV 71" 772 ^ 

l\ Q 2Q cot -S dx + 2j^ [ n \ Q 2Q cot -9 cos(2QA;x) = 

"2 + 2?? 2 fc=i-'-f+5zj 2 

1 + Sin ^j 1 tt 7J-J77, 

ln i ~iSr + n tan TT + V<3 2 ■ 

1 - sin ^ Q Q 

(26) 

Here the value of the last sum over k is estimated by integration by parts. It is 
of order 1/Q 2 . Thus, ( f2~4D reduces to 

r m As/Q 5 , s 

— = n (27) 

Qi 3 1 - cos ^ + sin ^ ' K J 

A similar analysis of (0) gives the result 

7T S I Tim 2 \ r m COS 2 ^ 

1 - tan — - — — — — — + 



2Qcos^ Q s \ Q 1-cos^ + sin^/ 2QP(1 + sin ^) ' 

(28) 
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From (^,^) we get 5 = ft = 1, s = n/2. The equation for the energy reads 
e = 4(1 - s/Q) + o(l/Q) = 4 - 27r/Q + o(l/Q). (Correspondingly, the lower 
boundary of the spectrum is e = —4 + 2ir/Q + o(l/Q).) The same correction 
2n/Q can be obtained from semi classical considerations |T8[ . 

4 Conclusions 

In this work we have studied the distribution of the roots Zk of the Bethe ansatz 
equations for the Harper equation in the limit as Q —>■ oo while the commen- 
surability parameter a — P/Q, P — 1. We considered two related sets of the 
Bethe ansatz equations to the main order as 1/Q — ► and obtained functions 
p(<f>) and x(0) that characterize the distribution. It is a continuous distribution 
of the roots on the unit circle. 

Considering the Bethe ansatz equations to the next (smaller) order in 1/Q 
we calculated the small a correction to the upper and lower boundaries of the 
spectrum of the Harper equation. 

One more step further could be to extend our treatment to the case of any 
fixed P, while Q —>■ oo. Numerical data suggests that in this case the roots Zk for 
the largest eigenvalue e max become continuously distributed on P smooth curves. 
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6 Appendix 

Here we derive the g-difference equation (||). 

Assume a = P/Q with coprime integers P, Q. Then (1) is the eigenvalue 
equation for a periodic tridiagonal matrix (denote this matrix M) with the period 
Q (that is the matrix elements = M i+ Qk+o). The approach to the spectral 
problem for such a matrix is well known [e.g., |I9 |. We look for solutions to ([3]) in 



the form ipj+Qk = e fij, j = 0, 1, . . . , Q — 1, where u is an arbitrary real number. 
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This reduces ([I]) to the following eigenvalue equation 
Mn-i + {q 2n e ie + q- 2n e~ ie )fi n + fi n+1 = efi n , 

n = 0, 1, . . . , Q - 1; q = e l7ra ; /i_i = e l > Q _i; /x Q = e*> - 

for a, Q x Q matrix (denote it L). It is easy to show that det(L — el) = 
(—1)Q(Sq(s;q>,9) — 2cosuj), where Sg(e;a,9) is a polynomial of degree Q in 
e, which is called the discriminant of M. The spectrum of M is obviously given 
by the equation Sq(e; a, 9) = 2 coscu, G [0, 7r]. As is known, the spectrum con- 
sists of exactly Q intervals, the image of [—2, 2] under the inverse of the transform 
Sq(e; a, 9) — A (see Figure |)|] Hence, all the extremal points of Sq(e; a, 6) are 
maxima and minima; moreover, at all maxima points Sg(:r max ; ot, 9) > 2 and at 
all minima points Sg(i min ; a, 9) < —2. 

As we mentioned in the introduction, the spectrum cr#(a) of the Azbel- 
Hofstadter problem is the union over all real 9 of the spectra of M. The Chambers 
formula 0] says that the polynomial Sg(e; a) = Sq(e; a, 9) + 2 cos 9Q does not 



depend on 9, which implies that cr^(a) is the image of [—4,4] under the inverse 
of the transform T,q(e; a) = A. Considering the Chambers formula at 9 = and 
9 = ir/Q, we see that EQ(a; max ; a) > 4 and Eg(x m i n ; a) < —4 at all its extremal 
points. So that cr^(a) also consists of Q intervals (see Figure £|). 

Let us now turn to the derivation of (0). Set uj = 0. Substituting /i n = 
Ejto U nk e ik6 <j) k , where f/ nfc = q 2nk , n, k = 0, 1, . . . , Q - 1,0 in we have 



fc _i + (g 2fe + g 2fe )0 fc + 4> k+1 = £4> k , 
fc = 0,l,...,g-l; 0_ x = e^Vo-H 
By the nondegenerate transformation 



(30) 



^=E? (fe+1/2)2+2fen Cn (31) 
n=0 



we get rid of the middle term in (|30|): 



(1 + q 2n )(n~l + (1 + g^ (n+1) )Cn + l = e( n , 

n = 0, 1, . . . , Q - 1; C-i = Co-i; Cq = Co- 



(32) 



6 One can prove this by showing that Sq(M; a, 8) — E, where the matrix elements of E are 
Eij = 1 if \i — j\ = Q and zero otherwise. The spectrum of E is obviously 2Q-degenerate, while 
the spectrum of the tridiagonal matrix M can be at most doubly degenerate. 

7 The matrix U is nondegenerate. det U ^ because UU* = QI. 
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Figure 4: The plot of Sq(e; a, 6) and Eq(s; a) for Q = 4 and some fixed 9. (For 
Q = 4 the curves do not depend on P: P = 1 or P = 3.) The thickest lines on the 
axis e are the intervals of the spectrum (7^(0:, 0) of M. The less thick lines are 
the intervals of the spectrum <7f/(a). Two of the latter intervals touch at e — 0. 
The curves shown are not the actual ones but qualitatively similar. 
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In order to obtain (|32| ) we had to assume the following condition: for Q even 
Q9 = 7r(mod 2n); for Q odd Q9 = 0(mod 2ir). If Q is even, we can make the 
shift ( n +Q/2 = Cm so that we get the equations 

(1 - Q 2n )U-i + (1 - q- 2{n+1) )Ui = 
n = 0,l,...,Q- 1 

with zero boundary conditions (£_i and £q may be arbitrary, e.g., zero). Further 
transformation £ n = p n i n q n ( n+1 ^ 2 gives 

*(? n - q~ n )Pn-i + ^(q n+1 - <T (b+1 W = ep n , 
n = 0,l,...,Q-l. 

Now multiplying ( |3~4"D by z n and then taking the sum over n from n = to 
Q — 1, we get 

iO" 1 + g;z)^(g,2) - iO" 1 + g^ 1 ^)^^^ 1 ^) = e^O), (35) 

which is equation @. Here ^(2) = Y^Zq z n p n - We derived ( ^5|) on condition 
that u = 0, Q6 = 7r(mod 2ir), and Q is even (hence, P is odd). It is easy to 
derive the same equation (|35|) for the two other possible cases (that is P even, 
Q odd, and P odd, Q odd) using slightly different from ([H|) transformations. In 
all cases we have the same restriction u — 0, Q# = 7r(mod 2tt). This implies that 
the solutions e of (0) are the points satisfying £q(£; a) = 0, so they lie inside the 
bands of the spectrum 07/ (a) of the Azbel-Hofstadter problem. As follows from 
the Chambers formula, they also lie in the bands of gm (ot, 9) for any 9. Take 
any band of ctm («,#)• It has a single point in it which is a solution of (Q). As 
9 changes from 9 = 0to9 = ir/Q the band moves through this point from one 
boundary where Sq(e; a, 0) = —2 to the other where Sq(e] a, ir/Q) = 2. 
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